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Abstract 

The chaotic synchronization between two bi-directionally coupled external 
cavity single-mode semiconductor lasers is investigated. Numerical simula- 
tion shows that anticipating synchronization and lag synchronization coexist 
in certain parameter regime. The anticipating time with different effects that 
were discussed quite differently in the previous theoretical analysis and exper- 
imental observation is determined by the involved parameters in the system. 
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I. Introduction 

In recent years, much attention has been paid to the chaotic synchronization in 
nonlinear systems [1-3]. There is much prospect of chaotic synchronization being 
applied in various ways, especially in secure communication where many works have 
been conducted and a lot of progress has been made [2,4,7]. In last few years, 
the theory of chaotic synchronization are being utilized in the research of neural 
networks [5]. 

Chaotic synchronization was introduced as studying the dynamics of several clas- 
sical models [1,2], such as Lorenz model. Then the experimental demonstration 
of the chaotic synchronization was given in the electrical systems and optical sys- 
tems successively [2,4,7]. Particular emphasis is put upon the synchronization in 
the chaotic external cavity semiconductor lasers [6-13] because of their ability to 
generate high-dimensional chaos and their ease of operation. Several years ago, 
the investigation in chaotic synchronization was confined to lag synchronizaiton [6]. 
There is a retardation time (r c ) between two lasers' output due to the finite time for 
the light to travel from the master to the slave. Recently, it is found that there is 
also anticipating synchronization [8-13] other than lag synchronization, the intensity 
of the slave is synchronized to the future intensity of the master. That means the 
slave can anticipate the dynamics of the master. 

Although anticipating synchronization is counterintuitive, its existence can be il- 
luminated both theoretically and experimentally [8-13]. Voss discovered anticipating 
synchronization in some simple models first [8]. Following this discovery, Masoller 
predicted anticipating synchronization in unidirectionally coupled lasers system [9] 
and showed that the anticipating time (ta) should be equal to the difference between 
round-trip time of the light in the master's external cavity and the retardation time, 
i.e., r a = r c — t. Recently Sivaprakasam reported the experimental demonstration 
of anticipating chaotic synchronization in a bi-directionally coupled external cavity 
semiconductor laser system [12] and found the anticipating time is irrespective of 
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the external-cavity round trip time r, i.e. r a = t c . At present there is no fully 
theoretical explanation of the difference between the theoretical anticipating time 
and the experimental observation [13]. 

In this paper, the synchronization between two bi-directionally coupled chaotic 
external cavity semiconductor lasers is investigated. Numerical simulations show 
very interesting results. (1) Anticipating synchronization and lag synchronization 
coexist in a system. Sometimes the system is inclined to anticipating synchroniza- 
tion, sometimes to lag synchronization. And long-term behavior is decided by the 
involved parameters such as the coupling strengths and the feedback rate in the 
master. (2) As an important character of the system, the anticipating time is also 
determined by the parameters. The parametric space is divided into several zones. 
In one of these zones, anticipating time is the same as derived in theoretic analysis 
reported in previous papers. In another zone, anticipating time is in agreement with 
experimental observation. 

II. Theoretical Model 

The dynamics of two bi-directionally coupled single-mode semiconductor lasers 
with only the master laser subjected to external optical feedback can be described 
by the widely utilized Lang-Kobayashi equations [14] 

dE n 
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Where subscript m and s denote the master and the slave respectively. E is the 
slowing varying complex field, and N is the normalized carrier density. The second 
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term on the right side in both (1) and (3) corresponds to the bi-directional coupling. 
k is the cavity loss, a the linewidth enhancement factor. Gi = iVj/(l + €i\Ei\ 2 ) is 
the optical gain where i stands for m or s and q is the gain saturation coefficient. 7 
is the feedback rate, r\ is the coupling strength and subscript sm indicates that the 
coupling is from the slave to the master, ms indicates that the coupling is from the 
master to the slave, ui is the optical frequency without feedback, Acu = u s — u m is 
the frequency detuning between the lasers, r is the external cavity round trip time 
in the master, and r c is the time for the light to fly from the master to the slave. £ 
is independent complex Gaussian white noise, and $ measures the noise intensity. 
j is the normalized injection current, and r is the carrier lifetime. For simplicity, 
here we take (3 m = (3 S = . 

III. Coexistence of lag synchronization and anticipating synchroniza- 
tion 

According to the equations(l)-(4), we make the numerical simulation of the dy- 
namics of two bi-directionally coupled semiconductor lasers. We figure out Q = 
cr 2 (arctan(j^)) to represent the synchronization quality, where I m and I s are the 
field intensities of the master and slave respectively (I m = \E m \ 2 ,I m = \E m \ 2 ), a 2 
is the variance calculation. Perfect synchronization is represented by Q = 0, on 
the other hand a high variance Q represents a poor synchronization. Because of 
the effect of anticipating synchronization or lag synchronization, we would derive a 
good synchronization represented by a low variance Q if the master laser output is 
shifted relatively to the slave output by an appropriate time t s . In this paper, t s 
takes t s > if the master laser output is shifted forward relative to the slave output. 
In this low variance Q represents an anticipating synchronization. t s takes 

t s < if the master laser output is shifted backward relative to the slave output, a 
low variance Q in this case represents a lag synchronization. 

We find that anticipating synchronization and lag synchronization coexist in the 
bi-directionally coupled semiconductor lasers system. The coexistence can be seen in 
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fig.l. Fig.l shows the dependence of the variance Q on the master laser time shift t s . 
there is a low variance on each side. The low variance on the right side(r s > 0) indi- 
cates the anticipating synchronization and the low variance on the left side(r s < 0) 
indicates the lag synchronization, the coexistence is further illustrated in detail in 
fig.2(a-l),2(a-2). Fig.2(a-1) and 2(a-2) plot the time traces of the field intensities of 
the master (upper) and the slave(lower). Each trace can be naturally divided into 
many small zones, labelled by mi, m 2 , m 3 , 777.4... in the master laser, si, s 2 , 83,54. ..in 
the slave laser, we detect the similarities between si and mi, m 2 and Si, s 2 and 
m 2 , m 3 and s 2 , s 3 and m 3 , m 4 and s 3 .... The dynamics of the system of two mu- 
tual coupled semiconductor lasers is of lag synchronization as far as the similarities 
between Si and m 1; s 2 and m 2 , s 3 and m 3 are concerned, on the other hand, the 
dynamics is of anticipating synchronization as far as the similarities between m 2 
and Si, m 3 and s 2 , m 4 and s 3 are concerned. Thus, the two kinds of synchroniza- 
tion coexist with each other. In addition, the similarities between mi and s«, s« 
and m i+ i fluctuate. The fluctuations indicate that sometimes the system behav- 
ior is inclined to anticipating synchronization orientation, sometimes is inclined to 
lag synchronization orientation. Six traces in fig. 2 share one dynamical trajectory. 
Fig.2(a-1) and 2(a-2) plot the intensities in the master and slave respectively in the 
range 250-300ns; fig.2(b-l) and 2(b-2) plot in the range 410-450ns, typically show- 
ing lag synchronization orientation; fig.2(c-l) and 2(c-2) plot in the range 900-930ns, 
typically showing anticipating synchronization orientation. 

Whether the long-term system behavior is anticipating or lag is decided by the 
involved parameters. In this paper the coupling strengths in both directions and the 
feedback rate in the master are discussed. We use AQ to indicate the orientation of 
the synchronization in a sufficient long term. 



Qa = 



min 



(T S > 0) 
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AQ = Q a - Qi 
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where a and 1 denote anticipating synchronization and lag synchronization respec- 
tively. The anticipating synchronization orientation is charactered by positive AQ 
(AQ > 0), while lag synchronization orientation is charactered by negative AQ 
(AQ < 0). 

In fig. 3 (a) we plot AQ as a function of the feedback rate in the master (^ m ) and 
the coupling strength from the slave to the master (fj sm ). The coupling strength from 
the master to the slave rj ms = 5.0. We can find AQ > in two regions labelled by ii, 
I2, where the long-term system behavior is inclined to anticipating synchronization 
orientation. In another region labelled by II, AQ < indicates the long-term system 
behavior is inclined to lag synchronization orientation. In region III, the long-term 
system dynamics is adiaphorous, represented by AQ 0. 

In fig. 3(b) we plot AQ as a function of the couplings in both directions(?7 ms ,^ sm ). 
We find that almost all the plane has AQ < 0, which indicates the long-term 
system behavior is inclined to lag synchronization in most cases. We mention that 
only the master has feedback which break the system's symmetry. The symmetry 
is considered in the following section. 

IV. Anticipating time 

There is a difference between the theoretical anticipating time reported in previ- 
ous papers and the experimental observation. Our numerical simulation shows that 
the two cases both occur in bi-directionally coupled chaotic external cavity semicon- 
ductor lasers system. In our numerical simulation, we take r c = 3.5ns, r = 6.7ns. 
Fig.3(c) plots the calculated anticipating time (r = t s ,t s > 0) as the function of 
the feedback rate in the master (^ m ) and the coupling strength from the slave to 
the master (f] sm ). In terms of r a one can clearly distinguish two regions. In the 
region labelled by I, r a is about 3.5ns, in agreement with experimental anticipating 
time (r a = r c = 3.5ns). In the region labelled by II in fig.3(c), r a is about 3.2ns, 
corresponding to previous theoretical reports (r a = r — r c = 3.2ns). 

Fig. 3(d) shows the dependence of anticipating time on the couplings in both 



6 



directions. The two regions, one for r a = r c = 3.5ns, the other for r a = r— r c = 3.2ns 
could also be distinguished clearly. 

With respect to the anticipating time, an important characteristic scalar, the 
system is symmetric. (?7 ms ,?7 sm ) = (a,b) and (77 ) = (b,a) give the same r a . The 
symmetry tell that the two lasers should be regarded as segments of a whole sys- 
tem. Two lasers participate in the dynamical evolution jointly, which could also be 
detected in fig. 2. We mention that the bi-directional coupling leads to the symmetry. 
And the symmetry indicates that the couplings in two directions are equivalent. 
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Figure Captions 



Fig. 1. Variance a 2 (arctan(j^-)) as a function of the shift time r s . The param- 
eters are k m = k s = 500ns, a m = a s = 3, e m = e s = 0.1, r nm = r ns = 1.0ns, 
jm = is = 101, cu m = io s = 0.3rad- ns -1 , 7 m = 2ns _1 , r^ ms = 5ns~\ r] sm = 3ns" 1 , 
r = 6.7ns, r c = 3.5ns 

Fig. 2. The time traces of the master laser field intensity (a-l)(b-l)(c-l) and 
the slave laser field intensity (a-2)(b-2)(c-2). (a-l)(a-2) range 250-300ns; (b-l)(b- 
2) range 410-450ns; (c-l)(c-2) range 900-930ns. 7 m = 10ns -1 , i] ms = 10ns -1 , 
i] sm = 2ns -1 [all other parameters are the same as fig.l]. 

Fig. 3. (a) The calculated AQ as a function of the feedback rate in the mas- 
ter 7 m and the coupling strength from the master to the slave r\ sm . i] sm = 5ns -1 . 
(b) The calculated AQ as a function of the couplings in both directions (r] ms , i] sm ). 
7m = 5ns -1 . (c) The calculated anticipating time as a function of 7 m and i] sm . 
Vsm = 5ns 1 . (d) The calculated anticipating time as a function of Tj ms and Tj sm . [all 
other parameters are the same as fig.l]. 
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